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ABSTRACT. Standard accretion disk theory is for- 
mulated which is based on the local heat balance. 
The energy produced by a turbulent viscous heating 
is supposed to be emitted to the sides of the disc. 
Sources of turbulence in the accretion disc are con- 
nected with nonlinear hydrodynamic instability, con- 
vection, and magnetic field. In standard theory there 
are two branches of solution, optically thick, and opti- 
cally thin. Advection in accretion disks is described by 
the differential equations what makes the theory nonlo- 
cal. Low-luminous optically thin accretion disc model 
with advection at some suggestions may become advec- 
tively dominated, carrying almost all the energy inside 
the black hole. The proper account of magnetic filed 
in the process of accretion limits the energy advected 
into a black hole, efficiency of accretion should exceed 
~ 1/4 of the standard accretion disk model efficiency. 
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1. Introduction 

Accretion is a main source of energy in binary X-ray 
sources, quasars and active galactic nuclei (AGN). The 
intensive development of accretion theory began after 
discovery of X-ray sources and quasars. Accretion into 
stars is ended by a collision with an inner boundary, 
which may be a stellar surface, or outer boundary of a 
magnetosphere for strongly magnetized stars. All grav- 
itational energy of the falling matter is transformed 
into heat and radiated outward. 

In black holes matter is falling to the horizon, from 
where no radiation arrives. All luminosity is formed on 
the way to it. The efficiency of accretion is not known 
from the beginning, and depends on angular momen- 
tum of the falling matter, and magnetic field embedded 
into it. It was first shown by Schwartsman (1971) that 
during spherical accretion of nonmagnetized gas the 
efficiency may be as small as 10^* for sufficiently low 
mass fluxes. He had shown that presence of magnetic 
field in the accreting matter increase the efficiency up 
to about 10%, and account of heating of matter due 
to magnetic field annihilation rises the efficiency up to 



about 30% (Bisnovatyi-Kogan and Ruzmaikin, 1974) 
In the case of a thin disc accretion, when matter has 
large angular momentum, the efficiency is about 1/2 
of the efficiency of accretion into a star with a radius 
equal to the radius of the last stable orbit. In the 
case of geometrically thick and optically thin accretion 
discs the situation is approaching the case of spherical 
symmetry, where presence of a magnetic field plays a 
critical role. 

Advection dominated accretion flow (ADAF) was 
suggested by Narayan and Yu (1995), and used as a 
solution for some astrophysical problems. The sugges- 
tions underlying ADAF: ignorance of the magnetic field 
annihilation in heating of a plasma fiow, electron heat- 
ing only due to binary collisions with protons (ions) 
had been critically analyzed in papers of Bisnovatyi- 
Kogan and Lovelace (1997, 1999), Bisnovatyi-Kogan 
(1999), Quataert (1997). There are contradictions be- 
tween ADAF model and observational data in radioe- 
mission of elliptical galaxies (Di Matteo et al., 1999), 
and X-ray emission of Seyfert galaxy NGC4258 (Can- 
nizzo, 1998). Account of processes connected with a 
small-scale magnetic field in accretion flow, strongly 
restricts solution. Namely, the efficiency of the accre- 
tion flow cannot become less then about 1/4 of the 
standard accretion disc value. 

2. Development of the standard model 

Matter falling into a black hole forms a disc when 
its angular momentum is sufficiently high. It happens 
when matter comes from the star companion in the 
binary, or from a tidal disruption of the star which tra- 
jectory approaches close to the black hole. The first 
situation is observed in many galactic X-ray sources 
(Cherepashchuk, 1996). A tidal disruption happens in 
quasars and active galactic nuclei (AGN) in the model 
of supermassive black hole surrounded by a dense stel- 
lar cluster. 

Equations of a standard accretion disk theory had 
been first formulated by (Shakura, 1972); some cor- 
rections and generalization to general relativity (GR) 
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had been done by Novikov and Thorne (1973). Ob- 
servational aspects of accretion disks have been ana- 
lyzed by Shakura and Sunyaev (1973). Note, that all 
authors of the accretion disc theory from USSR were 
students (N.I. Shakura) or collaborators (I. D. Novikov 
and R.A. Sunyaev) of academician Ya.B.Zeldovich, who 
was not among the authors, but whose influence on 
them hardly could be overestimated. The main idea of 
this theory is to describe a geometrically thin non-self- 
gravitating disc of a mass M^i, much smaller then the 
mass of a black hole M, by hydrodynamic equations 
averaged over the disc thickness 2h. 

2.1. Equations 

The small thickness of the disc in comparison with 
its radius h <^ r means small importance of the pres- 
sure gradient VP in comparison with gravity and in- 
ertia forces. Radial equilibrium equation in a disc is a 
balance between the last two forces with an angular ve- 
locity equals to the keplerian one = ^Ik = (^p^)^^^. 
Note, that just before a last stable orbit around a black 
hole this suggestion fails, but in the "standard" accre- 
tion disc model this relation is supposed to hold all 
over the disc, with an inner boundary at the last sta- 
ble orbit. The equilibrium equation in the vertical z- 
direction is determined by a balance between the gravi- 
tational force and pressure gradient ^ = —p '^^i^ . For 
a thin disc this differential equation is substituted by 
an algebraic one, determining the half-thickness of the 
disc in the form 
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The balance of angular momentum, related to the </> 
component of the Euler equation has an integral in a 
stationary case, written as 



M{j - jin) = -27rr^ 2htr^, t^ 



w-r- (2) 
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Here j = v^r = flr^ is a specific angular momen- 
tum, tr^ is a component of the viscous stress tensor, 
M > is a mass flux per unit time into a black hole, 
jin is equal to the specific angular momentum of mat- 
ter falling into a black hole. In the standard theory the 
value of jin is determined from physical considerations. 
For accretion into a black hole it is suggested, that on 
the last stable orbit the gradient of the angular veloc- 
ity is zero, corresponding to zero viscous momentum 
flux. In that case jin = O.Krf^, corresponding to the 
Keplerian angular momentum of the matter on the last 
stable orbit. 

The choice of the viscosity coefficient is the most 
speculative problem of the theory. In the laminar case 



at low microscopic (atomic or plasma) viscosity the sta- 
tionary accretion disc must be very massive and thick, 
and before its formation the matter is collected by disc 
leading to a small flux inside. It contradicts to obser- 
vations of X-ray binaries, where a considerable matter 
flux along the accretion disc may be explained only 
when viscosity coefficient is much larger. It was sug- 
gested by Shakura (1972), that matter in the disc is 
turbulent, what determines a turbulent viscous stress 
tensor, parametrized by a pressure 
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where Vs is a sound speed in the matter. This simple 
presentation comes out from a relation for a turbulent 
viscosity coefficient rjt ~ pvtl with an average turbulent 
velocity vt and mean free path of the turbulent element 
I. It follows from the definition of i,.^ in (^, when we 
take I « h from (n\) 



tr^ = pvthr 



dr 



PVtVs 



-apv. 
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where a coefficient a < 1 is connecting the turbulent 
and sound speeds vt = avg. Presentations of tr^ in 
(H) and (^ are equivalent, and only when the angu- 
lar velocity differs considerably from the Keplerian one 
the first relation to the right in (||) is more preferable. 
That does not appear in the standard theory, but may 
happen when advective terms are included. 

Development of a turbulence in the accretion disc 
cannot be justified simply, because a Keplerian disc 
is stable in linear approximation to the development 
of perturbations. It was suggested by Ya.B.Zeldovich, 
that in presence of very large Reynolds number Re = 
^ the amplitude of perturbations at which nonlin- 
ear effects become important is very low, so a turbu- 
lence may be developed due to nonlinear instability 
even when the disc is stable in linear approximation. 
Viscous stresses may arise from a magnetic field, it was 
suggested by (Shakura, 1972), that magnetic stresses 
cannot exceed the turbulent ones. It was shown by 
Bisnovatyi-Kogan and Blinnikov (1976), that inner re- 
gions of a highly luminous accretion discs where pres- 
sure is dominated by radiation, are unstable to vertical 
convection. Development of this convection produce a 
turbulence, needed for a high viscosity. 

With alpha- prescription of viscosity the equation of 
angular momentum conservation is written in the plane 
of a disc as 



Mij ~j,n)^'^Trr^aPoh. 



(5) 



When angular velocity is far from Keplerian one the 
relation (^) is valid with a coefficient of a turbulent 
viscosity 77 = apoVsoh, where values with the index 
"0" denote the plane of the disc. 

In the standard theory a heat balance is local, all 
heat produced by viscosity in the ring between r and 
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r + dr is radiated through the sides of disc at the same 
r. The heat production rate related to the surface 



unit of the disc is written as 



^ , dn 3 -GM 
Q+ = htr4,r— = —M^— 
dr OTT 



-j^ Jin 
j 



(6) 



Heat losses by a disc depend on its optical depth. The 
standard disc model (Shakura, 1972) considered a ge- 
ometrically thin disc as an optically thick in a vertical 
direction. That implies energy losses Q- from the disc 
due to a radiative conductivity, after a substitution of 
the differential equation of a heat transfer by an alge- 
braic relation 



Q- 



AacT* 
3 kS 
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Here a is a constant of a radiation energy density, c is a 
speed of light, T is a temperature in the disc plane, k is 
a matter opacity, and a surface density E — 2ph. Here 
and below p, T, P without the index "0" are related 
to the disc plane. The heat balance equation is repre- 
sented by a relation (5_|_ = Q^. A continuity equation 
in the standard model is used for finding of a radial 
velocity Vr 
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(8) 

Airrhp 2'KrYi 

Completing these equations by an equation of state 
P{p,T) and relation for the opacity k — k{p,T) we 
get a full set of equations for a standard disc model. 
For power low equations of state of an ideal gas P = 
Pg = pTZT {TZ is a gas constant), or radiation pressure 
P = Pr = and opacity in the form of electron 

scattering Kg, or Karammers formulae k^, the solution 
of a standard disc accretion theory is obtained analyti- 
cally. Checking the suggestion of a large optical thick- 
ness confirms a self-consistency of the model. Note 
that solutions for different regions of the disc with dif- 
ferent equation of states and opacities are not matched 
to each other. 



2.2. Optically thin solution 

It was found by Shapiro ct al. (1976) that there 
is another branch of the solution for a disc structure 
with the same input parameters M, M, a which is also 
self-consistent but has a small optical thickness. That 
implies another equation of energy losses, determined 
by a volume emission w qph. The emissivity of 
the unit of a volume q is connected with a Planckian 
averaged opacity Kp by a relation q « acT^Kp. In the 
optically thin limit the pressure is determined by a gas 
P — Pg. In the optically thin solution the thickness of 
the disc is larger then in the optically thick one, and 
density is lower. 



While heating by viscosity is determined mainly by 
heavy ions, and cooling is determined by electrons, the 
rate of the energy exchange between them is important 
for a structure of the disc. The energy balance equa- 
tions are written separately for ions and electrons. For 
small accretion rates and lower matter density the rate 
of energy exchange due to binary collisions is so slow, 
that in the thermal balance the ions are much hot- 
ter then the electrons. That also implies a high disc 
thickness. In the highly turbulent plasma the energy 
exchange between ions and electrons may be strongly 
enhanced due to presence of fluctuating electrical fields, 
where electrons and ions gain the same energy. In such 
conditions difference of temperatures between ions and 
electrons may be negligible. The theory of relaxation 
in the turbulent plasma is not completed, but there are 
indications to a large enhancement of the relaxation in 
presence of plasma turbulence, in comparison with the 
binary collisions (Galeev and Sagdeev, 1983; Quataert, 
1997). 

2.S. Accretion disc structure from equations describ- 
ing continuously optically thin and thick regions 

Equations of the disc structure smoothly describing 
transition between optically thick and optically thin 
disc, had been obtained using Eddington approxima- 
tion. The expressions had been obtained (Artemova et 
al., 1996) for the vertical energy flux from the disk Fq, 
and radiation pressure in the symmetry plane 
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, • At To > > 1 we have (0) from @. 
In the optically thin limit ^ tq <C 1 we get 



Fq — acT^TaO, Prad,{ 



-acT^Tao. 



(10) 



Using Fq instead of and equation of state P = 
pTZT + Prad,o, the equations of accretion disc struc- 
ture together with equation — Fq, with from 
(|^), have been solved numerically by Artemova et al. 
(1996). Two solutions, optically thick and thin, exist 
separately when luminosity is not very large. They in- 
tersect at TO = TOf, and there is no global solution for ac- 
cretion disc at TO > TOf,. It was concluded by Artemova 
et al. (1996) that in order to obtain a global physically 
meaningful solution at to > ifib, account of advectiv- 
ion is needed. For the calculated case Mbh = 10^ Mq, 
a = 1.0 at TO = mi, luminosity of the accretion disk is 
less than the critical Eddington one. 

3. Accretion discs with advection 
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Standard model gives somewhat nonphysical behav- 
ior near the inner edge of the accretion disc around 
a black hole, with a zero heat production at the in- 
ner edge of the disk. It is clear from physical ground, 
that in this case the heat brought by radial motion of 
matter along the accretion disc becomes important. In 
presence of this advective heating (or cooling) term, 
depending on the radial entropy S gradient, written 
as Qadv = ^^^) the equation of a heat balance is 
modified to 



Q+ + Qadv — Q- 



(11) 



In order to describe self-consistently the structure of 
the accretion disc we should also modify the radial disc 
equilibrium, including pressure and inertia terms 



1 dP 



dVr 



p dr 



dr 



(12) 



Appearance of inertia term leads to transonic radial 
flow with a singular point. Conditions of a contin- 
uous passing of the solution through a critical point 
choose a unique value of the integration constant jin- 
First approximate solution for the advective disc struc- 
ture have been obtained by Paczyhski and Bisnovatyi- 
Kogan (1981). Attempts to find a solution for advec- 
tive disc had been done by Abramovicz et al. (1988), 
Matsumoto et al. (1984). For moderate values of 
M a unique continuous transonic solution was found, 
passing through singular points, and corresponding to 
a unique value of jin. The number of critical point 
in the radial flow with the gravitational potential 0g 
(Paczynski and Wiita, 1980) (^g = Vg = 

may exceed unity. Appearance of two critical points 
for a radial flow in this potential was analyzed by 
Chakrabarti and Molteni (1993). Using of equations 
averaged over a thickness of the disc changes a struc- 
ture of hydrodynamic equations, leading to a position 
of singular points not coinciding with a unit Mach num- 
ber point. 



Kogan and Ruzmaikin (1974). For a spherical accre- 
tion with V = {vr, 0, 0) the equations describing a fleld 
amplification in the ideally conducting plasma reduce 
to (Bisnovatyi-Kogan, Ruzmaikin, 1974) 



d{r^Br) 
di 



= 0, 



djrVrBe) 
di 



= 0, 



d{rVrB^) 
di 



= 0, (13) 



where ^ = ■§t + v^-^ is a full Lagrangian deriva- 



tive. Consider a free fall 



with 



2GM 



The initial condition problem is solved separately for 
poloidal and toroidal fields. For initially uniform field 
Bro = BqcosO. BgQ = — Sosinf? we get the solution 
(Bisnovatyi-Kogan, Ruzmaikin, 1974) 



= Booo^O^.rs^ ^ _:^$V3, (14) 



where $i — r^/^ + |t\/2GM. The radial component of 
the field is growing most rapidly. It is ~ for large 
times, ~ 

^4/3 

at given small radius, and is growing with 
time everywhere. For initially dipole magnetic field 



BrO — 



Bq cos 6 



90 — 



Bo sin 6 
2r3 



we obtain the following solution 
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, Be 
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(15) 



Here the magnetic field is decreasing everywhere with 
time, tending to zero. That describes a pressing of 
a dipole magnetic field to a stellar surface. The az- 
imuthal stellar magnetic field if confined inside the 
star. When outer layers of the star are compressing 
with a free-fall speed, then for initial field distribution 
B^o = Bor"- f{9) the change of B^ with time is de- 
scribed by a relation B^ = — 

5. Two-temperature advective discs 



4. Amplification of the magnetic field at a 
spherical accretion 

A matter flowing into a black hole is usually mag- 
netized. Due to more rapid increase of a magnetic en- 
ergy the role of the magnetic field increases when mat- 
ter flows inside. It was shown by Schwartsman (1971), 
that magnetic energy density Em approaches a density 
of a kinetic energy E}~, and he proposed a hypothesis 
of equipartition Em ~ Ek, supported by continuous 
annihilation of the magnetic fleld in a region of the 
main energy production. This hypothesis is usually 
accepted in the modern picture of accretion (Narayan 
and Yu, 1995). Account of the heating of matter by 
magnetic fleld annihilation was done by Bisnovatyi- 



In the optically thin accretion discs at low mass 
fluxes the density of the matter is low and energy ex- 
change between electrons and ions due to binary colli- 
sions is slow. In this situation, due to different mech- 
anisms of heating and cooling for electrons and ions, 
they may have different temperatures. First it was real- 
ized by Shapiro et al. (1976), where advection was not 
included. It was noticed by Narayan and Yu (1995), 
that advection in this case is becoming extremely im- 
portant. It may carry the main energy flux into a black 
hole, leaving rather low efficiency of the accretion up 
to 10^'^ — 10~^ (advective dominated accretion flows 
- ADAF). This conclusion is valid only when the ef- 
fects, connected with heating of matter by magnetic 
fleld annihilation are neglected. 
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In the ADAF solution the ion temperature is about 
a virial one kTi ~ GMmi/r, what means that even 
at high initial angular momentum the disc becomes 
thick, forming a quasi-spherical accretion flow. When 
energy losses by ions are low, some kind of a "thermo- 
viscous" instability is developed, because heating in- 
creases a viscosity, and viscosity increases a heating. 
Development of this instability leads to formation of 
ADAF. 

A full account of the processes, connected with a 
presence of magnetic field in the flow, is changing 
considerably the picture of ADAF. It was shown by 
Schwarzman (1971), that in the region of the main en- 
ergy production, the condition of equipartition takes 
place, and efficiency of a radiation increase enor- 
mously from ^ 10~* up to ~ 0.1 due to magneto- 
bremstrahlung. To support the condition of equiparti- 
tion a continuous magnetic field reconnection and heat- 
ing of matter due to Ohmic dissipation takes place. 
It was shown by Bisnovatyi-Kogan and Ruzmaikin 
(1974), that due to Ohmic heating the efficiency of a 
radial accretion into a black hole may become as high 
as ^ 30%. The rate of the Ohmic heating in the con- 
dition of equipartition was obtained in the form 



dr 



3 



2 87rpr 



(16) 



In the supersonic flow of the radial accretion equipar- 
tition was suggested in a form (Schwartsman, 1971) 

^ Ki ^ — P^^' . For the disc accretion equiparti- 
tion between magnetic and turbulent energy was sug- 
gested by Shakura (1972), what reduces with account 
of " alpha" prescription of viscosity to a relation ^ ^ 

2 

^ « Q!^P, where am characterizes a magnetic viscos- 
ity in a way similar to the turbulent a viscosity. It was 
suggested by Bisnovatyi-Kogan and Ruzmaikin (1976) 
the similarity between viscous and magnetic Reynolds 
numbers, or between turbulent and magnetic viscosity 
coefficients Re = RCm = where the turbu- 

lent magnetic viscosity rjm is connected with a turbu- 

2 

lent conductivity a ~ . Taking rim ~ ^^Vj get 
a turbulent conductivity 



(17) 



in the optically thin discs. For the radial accretion the 
turbulent conductivity may contain mean free path of 
a turbulent element k in ( |l7|) instead of h. In ADAF 
solutions, where ionic temperature is of the order of the 
virial one two above suggestions for magnetic equipar- 
tition almost coincide at am ~ 1- 

The heating of the matter due to an Ohmic dissi- 
pation may be obtained from the Ohm's law in radial 
accretion 



T— = « -o— 



dr 



pVr 



B^vl 



pVrC^ 



AnpamVrVsk 



, (18) 



what coincides with ( |l^ ) when a™ — 3^-^77^, or /j = 
„ ^''"-^ — . Here a local electrical field strength in a 

highly conducting plasma is of the order of £ ~ ^^7^1 
ve ^ vt ^ avs for a radial accretion. 

Equations for a radial temperature dependence in 
the accretion disc, separately for the ions and electrons 
are written as 

— rr ^~J7 ~ '^ni + T^Bi — Qie , (19) 

dt dt 



dEf, Pf, dp 



'Hse + Qie — Cff — Ccyc , (20) 



dt p2 

A rate of a viscous heating of ions 7i,,i is obtained from 
i) as 



27rr 



3 vkv; 

= 7:Oi 

2 r 



T^rie < -1/ Ti-rji- (21) 



Combining (1) , (|) , (|) , we get 



a 



vkJ' 



V2^r, 

VK 



M vlJ 
47ra\/2 r'^vj 



(22) 

where vk = rflx, J = 1 ~ The rate of the 

energy exchange between ions and electrons due to 
the binary collisions was obtained by Landau (1937). 
Neglecting pair formation for a low density accretion 
disc, we may write an exact expression for a pres- 
sure Pg = P, + Pe ^ n^kTe + UekTp = nek{Te + Ti), 

and an approximate expression for an energy, con- 
taining a smooth interpolation between nonrelativis- 
tic and relativistic electrons. The bremstrahlung C// 
and magneto-bremstrahlung Ccyc cooling of maxwellian 
semi-relativistic electrons, with account of free-bound 
radiation in nonrelativistic limit, may be written as 
interpolation of limiting cases (Bisnovatyi-Kogan and 
Ruzmaikin, 1976). 

In the case of a disk accretion there are several char- 
acteristic velocities, vk, Vr, Vg, and Vt = avg, all of 
which may be used for determining "equipartition" 
magnetic energy, and one characteristic length h. Con- 
sider three possible choices of w^=w|-, w^, and for 
scaling B^ = Airpvg. The expression for an Ohmic 
heating in the turbulent accretion disc also may be 
written in different ways, using different velocities ve in 
the expression for an effective electrical field £ — 
A self-consistency of the model requires, that expres- 
sions for a magnetic heating of the matter Ti.B, ob- 
tained from the condition of stationarity of the flow 
(p^), and from the Ohm's law (^8|), should be iden- 
tical. It happens at j^-^ — That imphes 
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ve ^ Vr ^ °''"%T — '"^^hr < Vt- In the advec- 
tive models J is substituted by a function which is 
not zero at the inner edge of the disc. The heating 
due to magnetic field reconnection Ti^ in the equa- 
tions (|l^) , ( pO| ) may be written with account of (^l|) as 

= 4^^vr = ^H^, (^)'. At VB = VK the ex- 
pressions for viscous and magnetic heating are almost 
identical. Observations of the magnetic field recon- 
nection in the solar flares show (Tsuneta, 1996), that 
electronic heating prevails over the ionc one. Transfor- 
mation of the magnetic energy into a heat is connected 
with the change of the magnetic flux, generation of the 
vortex electrical field, accelerating the particles. 

The equations (19), ( |2C| ) have been solved by 
Bisnovatyi-Kogan and Lovelace (1997) for nonrela- 
tivistic electrons, at vb—vk- The combined heat- 
ing of the electrons and ions were taken as Tie ~ 
(2 — g)'Hni, We = g'Hni- The results of calculations 
for g ~ 0.5 1 show that almost all energy of the 
electrons is radiated, so the relative efficiency of the 
two-temperature, optically thin disc accretion cannot 
become lower then 0.25. Increase of the term Qif, due 
to plasma turbulence may restore the relative efficiency 
to a value, corresponding to the optically thick disc. 



6. Discussion 

Observational evidences for existence of black holes 
inside our Galaxy and in the active galactic nuclei 
(Cherepashchuk, 1996; Ho, 1999) make necessary to 
revise theoretical models of the disc accretion. The 
improvements of a model are connected with account 
of advective terms and more accurate treatment of the 
magnetic field effects. Account of the effects connected 
with magnetic field annihilation does not permit to 
make a relative efficiency of the accretion lower then 
~ 0.25 from the standard value. Strong relaxation con- 
nected with the plasma turbulence may increase the 
efficiency, making it close to unity. For explanation of 
underluminous galactic nuclei two possible ways may 
be suggested. One is based on a more accurate esti- 
mations of the accretion mass flow into the black hole, 
which could be overestimated. The second is based on 
existence of another mechanisms of the energy losses 
in the form of accelerated particles, like in the radio- 
pulsars, where these losses exceed strongly a radiation. 
This is very probable to happen in a presence of a 
large scale magnetic field which may be also respon- 
sible for a formation of the observed jets (Bisnovatyi- 
Kogan, 1999; Blandford and Begelman, 1999). We may 
suggest, that underlumilnous AGN loose main part of 
their energy to the formation of jets, like in SS 433. 
The search of the correlation between existence of jets 
and lack of the luminosity could be very informative. 
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